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1. Introduction 

There is much interest in precision lattice calculations of K — > %% decays because they can 
yield information on the origin of the AI = I rule and CP violation in the Standard Model [jl|, 
While many quenched calculations of these decays have been performed (for example [Q, ||, ^ |5|, 
|6]]), a full dynamical calculation has yet to be done. Steps have been made here toward a complete 
calculation of this type. In order to obtain reasonable precision we use 2+1 flavors of domain wall 
fermions (DWF) on a 24 3 x 64 or even 32 3 x 64, L s = 16 lattice. 

In physical K — > %% decays, the kinematics are such that the pions have non-zero momentum 
in the CM frame. Giving the pions momentum on the lattice can introduce a lot of noise due to 
the fact that one is projecting onto an excited state of the two pion system rather than the ground 
state. For this reason only pions with (nearly) zero momentum are simulated in this paper, and the 
introduction of significant pion momentum is left to future work. 

2. Four Quark Operators and the Effective Hamiltonian 

The weak interactions and the effects of the heavier quarks can be included in the lattice QCD 
simulation by evaluating matrix elements of an effective Hamiltonian [|7[ ||]. In particular we use 
the conventions of [0J 

Q 10 

^ AS=1 = -^V^V*£[z / (M) + Ty 1 -(M)]G« (2-1) 
V2 I=1 

where Vu are CKM matrix elements, Zi and y,- are Wilson coefficients, X = — A,/A„ with Xj = Vj c iVJ s , 
and {Qi,i = 1, 10} are four quark operators. Therefore we are interested in calculating matrix 
elements of the four quark operators Qj between a K and a %% state. These operators can be split 
into AI = 3/2 and AI = 1/2 parts, where AI is the change in isospin induced by the operator. They 
can then be further classified by how they transform under the chiral SU (3)l x SU (3)r symmetry, 
and the representations (27,1), (8,8), and (8,1) are all found among various of the operators [|l|, 0]. 

3. Extraction of the Matrix Element on the Lattice 

We calculate only matrix elements of the operator Q^jlj 2 , the single operator that transforms 
as (27,1) under SU(3)l x SU(3)r and AI = 3/2 under isospin. To simplify matters we calculate 
the unphysical matrix element (tc + tc + \Q'^!,~^ 2 \K + ) which can be related to the physical matrix 

element {n + n°\Q^J 2 \K+) by the Wi gner Eckhart theorem ^ if Q'^j i) i s given by 

e'Sf = s '^( l - - f)d (3.i) 

To extract the matrix element of this operator we calculate the following correlation functions 

C K (t,t K ) = (0 K (t)O f K (t K )), C nn (t,t n ) = (0jtw(O0Li(ljO), Cff(t n ,t,t K ) = {O njt (t„)0 w (t)Ol(t K )) 

(3.2) 

where tx <t <t n and where the interpolating operators are given by 

Oi(0 = E«(x,0r 5 *(y»0 J <W0= L d(x,t)y s u(y,t)d(x',t)r i u(y',t) (3.3) 

x,y x.y.x'^' 

w (t) = £j(x,07m(1 - 7 5 M(x,/)«(x,0r /i (l - f)d(x,t) (3.4) 
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Note that the interpolating operators in (3.3) are such that the correlators use wall sources and wall 
sinks, and project onto the zero momentum kaon and nearly zero momentum two pion states. 
For tx -C t <C t n we expect the following quotient of correlators to show a plateau in t: 

t„,t,t K ) Jt t ...... 

t K <.t <t n (3.5) 



CK{t,tK)C7C7i{t,t7c) ZtckZk' 

Here M is the matrix element we wish to extract. Zk and Z %% appear in the normalization factors 
for the kaon and two pion correlators respectively, and it is only possible to extract 



4. Details of the Lattice Calculation 

Calculations were carried out on the RBC/UKQCD 32 3 x 64, L s = 16 2+1 flavor domain wall 
fermion lattices. The inverse lattice spacing for these lattices is 2.42(4) GeV corresponding to 
a physical volume of (2.6 fm) 3 . The sea strange quark mass is always 0.03 in lattice units, and for 
the sea light quark mass (henceforth denoted by just m sea ) there is an ensemble with m sea = 0.004 
and an ensemble with m sea = 0.008, both with 129 configurations. For each ensemble, inversions 
are performed with the following valence masses: m va i =0.002, 0.004, 0.006, 0.008, 0.025, and 
0.030. We calculate matrix elements for all possible valence mass combinations such that m s > mi 
where m s is the valence strange quark mass and mi is the valence light quark mass. 

We add and subtract propagators with periodic and antiperiodic boundary condtions in order 
to double the effective time length and suppress around the world contributions. The resultant 
periodic plus antiperiodic (P+A) propagator has a source at t=0 and the resultant periodic minus 
antiperiodic (P-A) propagator effectively has a source at t=64. These provide the left and right 
walls for the kaon at tg = and the two pions at t n = 64 respectively, and the time t at which the 
operator is located is varied. Twelve wall source propagators using a unit source distributed over a 
single time slice, each with a specific spin and color (fixed in Coulomb gauge), are computed on 



each configuration and used to evaluate the correlation functions in (3.2). 



5. Results 

Effective mass plots of the kaon and two pion correlators are shown in Figure [I], for m sea = 
0.004, m s = 0.03, and mi = 0.004. When the K —> nn correlator is divided by the kaon and two 



pion correlators as in (p.5|), we find that the quotient shows a plateau as expected. The quotient and 
a fit to the plateau are shown for m sea = 0.004, m s = 0.03, and m; = 0.004 in Figure ||. 

From the lattice matrix element we wish to find the physical quantity \A2\, as defined in 
[^], which can be compared to experiment. |./#| must be multiplied by the following factors in 
order to obtain |A2| 



1 . The Wilson coefficient plus factors that appear in the effective Hamiltonian (gjj) multiplying 
the operator. Wilson coefficients evaluated at pL = 2 GeV are interpolated from [JToj] . 

2. Numerical factors due to the Wigner Eckhart transformation that relates K + — ► n + 7l° to 

K+ -> 7T+7T+. 

3. A lattice to MS renormalization factor Z^ s (/x) for the operator. It is evaluated at the same 



scale /I as the Wilson coefficient and is obtained from [11] 
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Figure 1: Effective mass plots and fitted values of the mass/energy for the kaon and two pion correlators 
with m sea = 0.004, m s = 0.03, and nti — 0.004. Left: Kaon correlator. Right: Two pion correlator. 
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Figure 2: Plot of the K — ► Kit correlator divided by the kaon and two pion correlators for m sea = 0.004, 
m s = 0.03, and m; = 0.004. The fit to the plateau is shown. 



4. A factor proportional to y 'm^E'^JJ' ', where L = 32c? is the spatial extent of the lattice. This 
performs a naive transformation between finite and infinite volume normalization, ignoring 



the effects of %% scattering that have been analyzed in [12]. We hope to include the more 



accurate correction of [12] in future work 



Technically the contribution of the other two AI = 3/2 four quark operators to \A2\ should also 
be included. However, these have been found to be much smaller than the contribution of the (27, 1) 
operator [^], so we neglect them. 

A table of results for all of the different mass combinations can be found in Table 1. The 
data for [A2I is new and was not presented at the time of the talk. Note that none of the mass 
combinations correspond to a process in which energy is conserved, although 
0.03, mi = 0.004 comes close. Notice also that changing m sea doesn't have much of an effect. This 
insensitivity to the sea quark mass justifies doing a linear interpolation/extrapolation to the unitary 
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point m sea = mi with respect to the sea quark mass for a given ra/. The sea strange quark mass, 
however, is always fixed at 0.03, so we are stuck with a non-unitary sttange quark if m s ^ 0.03. 

Since some of the mass combinations come close to conserving energy, it is not unreasonable 
to attempt to interpolate/extrapolate to energy conserving kinematics. Specifically, holding the pion 
mass constant we plot [A2I as a function of m\ — E\ % {rtiK is varied, E n% is constant), fit to a straight 
line, and extrapolate or interpolate to m\ — E\ % = 0. We do this for four different pion masses, and 
the plot plus extrapolation for m n = 308 MeV is shown in Figure [| on the left. 

We now have values of \At\ for several pion masses extrapolated to kinematics which are 
energy conserving, as tabulated in Table 2. In Figure |3] on the right, we add the preliminary data in 
this work to the plot of |A 2 | ~ Re(A 2 ) vs. m\ found in [§]. 

In Figure ^| we see that the data in the present work does not appear to agree with that in [0]. 
This may be due to the fact that dynamical quarks are used in the present work whereas [||] was 
done in the quenched approximation. However, it could easily be due to the different methods used 
to set the lattice spacing (the p mass for the lattices used in ||6|] vs. the Q. mass for the lattices used 
in this work, see [^]), especially considering that IA2I is a dimensioned quantity that is proportional 
to a~ 3 making it extremely sensitive to errors in a. Furthermore, the lattice spacings in the two 
works are quite different (a -1 = 2.42 GeV in this work compared to aT x = 1.31 GeV in ^) so that 
finite lattice spacing effects could come into play. It has been checked that the same conventions 
for the normalization of IA2I were used in this work and in [0]. 
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Figure 3: Left: Plot of the physical quantity |A 2 | vs. m 



Pion Mass Squared (GeV ) 



E\ K for the fixed pion mass m K = 308 MeV, and 



fixed E nK — 626 MeV. This plot is fit with a straight line, and extrapolated to the energy conserving point 
m\ — E\ % = (red circle). Right: Plot of |A2 1 vs. m\. The preliminary data points of this work (red, dashed 
line) are shown alongside those from (black, solid line). For the two largest pion masses in the present 
work the extrapolations to energy conserving kinematics are very large. Thus the error bars for these masses 
are the systematic rather than statistical error, which is estimated as the difference between the extrapolated 
value obtained from a linear and from a quadratic fit. 



6. Conclusion 

We have calculated the contribution of the operator to 1^2 1 in A/ = 3/2 K — > 7171 

decays with momentum pions on a 32 3 x 64, L 4 = 16 2+1 flavor DWF lattice. We did a linear 
extrapolation to a unitary light sea quark mass from the two light sea quark masses of 0.004 and 
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Table 1: Kaon masses, two pion energies, matrix element values, and IA2I results for different valence mass 
combinations. Each line of data for m sea = 0.004 has the corresponding line of data for m sea = 0.008 below 
it so that the insensitivity of the data to sea quark mass can be seen. 



ffll 


m s 




mj( (MeV) 


Enir (MeV) 


\J(\ 


\A 2 \ (10~ 8 GeV) 


0.002 


0.002 


0.004 


237.1(9) 


485(2) 


0.001 166(27) 


1.510(36) 


0.002 


0.002 


0.008 


241.4(9) 


493(3) 


001 148(29) 


1.521(41) 


0.002 


0.004 


0.004 


274.9(9) 


485(2) 


001 189(28) 


1 657(40) 


0.002 


0.004 


0.008 


279 0(8) 

— is . \j \yj j 


493(3) 


001 159(28) 


1 651(42) 


0.002 


0.006 


0.004 


308 1(9) 


485(2) 


001213(30) 


1 790(45) 


0.002 


0.006 


0.008 


312.1(8) 


493(3) 


001 178(28) 


1.774(45) 


0.002 


0.008 


0.004 


338 1 (9) 


485(2) 


001236(31) 


1 910(50) 


0.002 


0.008 


0.008 


341.9(8) 


493(3) 


001200(29) 

v./ ■ v./ v./ 1 ^— v./ v./ y z— y / 


1 893(48) 


0.002 


0.025 


0.004 


528(1) 


485(2) 


001399(50) 


2 701(97) 


0.002 


0.025 


0.008 


531(1) 


493(3) 


001409(47) 

\j - \J \J i- vyj s \ r / / 


2 770(94) 


0.002 


0.03 


0.004 


572(1) 


485(2) 


001438(55) 


2.89(11) 


0.002 


0.03 


0.008 


576(1) 


493(3) 


001472(55) 

V.UU 1 1 / ^- V -y -J J 


3.01(12) 


0.004 


0.004 


0.004 


307.7(9) 


626(2) 


0.001459(29) 


2 773(56) 


0.004 


0.004 


0.008 


311.7(8) 


632(2) 


001408(26) 

v./ . v./ v./ 1 rv./ yj \ z_» v./ / 


2 721(51) 


0.004 


0.006 


0.004 


337 5(8) 


626(2) 

yr — yr \ — j 


001466(30) 


2 918(60) 

^ « .y 1 y yj y_/ j 


0.004 


0.006 


0.008 


341 3(8) 


632(2) 


001418(25) 


2 869(53) 

i- . U U ' y -y ) 


0.004 


0.008 


0.004 


364 8(8) 


626(2) 


001477(30) 


3 056(64) 


0.004 


0.008 


0.008 


368 6(8) 


632(2) 


001433(26) 


3 012(55) 


0.004 


0.025 


0.004 


545(1) 


626(2) 

yr — y*\ — / 


001609(41) 

yj *yjyj x. \y \y s \ r j. / 


4.07(10) 


0.004 


0.025 


0.008 


549.0(9) 


632(2) 


001 576(33) 

U.UUl -y 1 \J \ ^y ^y J 


4 044(87) 


0.004 


0.03 


0.004 


588(1) 


626(2) 


001 644(44) 


4.32(12) 


0.004 


0.03 


0.008 


592(1) 


632(2) 


001619(37) 


4.31(10) 


0.006 


0.006 


0.004 


364 6(8) 


739(2) 


001668(31) 


4 076(76) 

r m \y / \y \ 1 vy / 


0.006 


0.006 


0.008 


368.4(7) 


746(2) 


001623(26) 


4.023(64) 


0.006 


0.008 


0.004 


389 9(8) 

^y yj J . J \ yj J 


739(2) 


001673(31) 


4 228(79) 


0.006 


0.008 


0.008 


393.7(7) 


746(2) 


001633(26) 


4 182(66) 


0.006 


0.025 


0.004 


562 5(9) 


739(2) 


001767(37) 


5.37(11) 


0.006 


0.025 


0.008 


566 1(8) 

~y \j yj ■ a y yj § 


746(2) 


001746(31) 


5 362(95) 


0.006 


0.03 


0.004 


604 6(9) 


739(2) 

/ oy\y*} 


001798(39) 

U.UU 1 ( y 1 yJ\ y J 


5 66(12) 


0.006 


0.03 


0.008 


608 3(9) 


746(2) 


001780(33) 


5.67(11) 


0.008 


0.008 


0.004 


413.6(8) 


837(2) 


0.001863(32) 


5.492(95) 


0.008 


0.008 


0.008 


417.4(7) 


843(2) 


0.001815(26) 


5.414(80) 


0.008 


0.025 


0.004 


579.4(8) 


837(2) 


0.001937(35) 


6.76(13) 


0.008 


0.025 


0.008 


582.9(8) 


843(2) 


0.001900(30) 


6.70(11) 


0.008 


0.03 


0.004 


620.5(8) 


837(2) 


0.001964(37) 


7.09(14) 


0.008 


0.03 


0.008 


624.0(8) 


843(2) 


0.001928(32) 


7.03(12) 


0.025 


0.025 


0.004 


710.2(6) 


1427(1) 


0.003172(48) 


20.89(32) 


0.025 


0.025 


0.008 


712.6(7) 


1431(2) 


0.003101(42) 


20.51(28) 


0.025 


0.03 


0.004 


745.1(6) 


1427(1) 


0.003180(48) 


21.45(33) 


0.025 


0.03 


0.008 


747.3(7) 


1431(2) 


0.003115(43) 


21.09(29) 


0.03 


0.03 


0.004 


778.7(6) 


1564(1) 


0.003509(53) 


26.52(40) 


0.03 


0.03 


0.008 


780.6(7) 


1567(8) 


0.00344(13) 


26.1(1.1) 
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Table 2: \A%\ vs. m\ after an interpolation/extrapolation to energy conserving kinematics. 



ml (GeV 2 ) 


|A 2 | (10~ 8 GeV) 


0.05521(68) 
0.09470(52) 
0.13432(40) 
0.17425(61) 


2.46(14) 
4.64(14) 
6.96(12) 
9.50(20) 



0.008 studied in our calculation, but did not attempt to make the sea and valence strange quark 
masses agree. We interpolated to energy conserving kinematics as well. When this preliminary data 
is plotted alongside the data of ^ in a graph of |A2| vs. m 2 n we see a disagreement which could be 
explained by different methods of determining the lattice spacing and finite lattice spacing errors. 
Also, the present work has the advantage of including 2+1 flavors of dynamical quarks, whereas 
all previous work has been in the quenched approximation. In the present work there is only data 
for decays to zero momentum pions. 

Future plans include doing the full AI = 3/2 calculation on 32 3 x 64 lattices, but with signifi- 
cantly stronger coupling so that the resulting larger physical volume will permit using pion masses 
much closer to the physical value. Momentum will be given to the pions using twisted boundary 
conditions [|3| |], ||, 13] and the kinematics will be very close to physical kinematics. 
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